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I ■ ABSTRACT 

> : 

, Context. Young massive clusters are usually deeply embedded in dust and gas. They represent perfect astrophysical laboratories 
■ for study of massive stars. Clusters with Wolf-Rayet (WR) stars are of special importance, since this enables us to study a 
coeval WR population at a uniform metallicity and known age. 

Aims. We started a long-term project to search the inner Milky Way for hidden star clusters and to study them in details. 
00 GLIMPSE 30 (G30) is one of them. The cluster is situated near the Galactic plane (;=298?756, 6=-0?408) and we aimed to 
determine its physical parameters and to investigate its high-mass stellar content and especially WR stars. 
1 Methods. Our analysis is based on SOFI/NTT JsHKs imaging and low resolution (R~2000) spectroscopy of the brightest 
' cluster members in the K atmospheric window. For the age determination we applied isochrone fits for MS and Pre-MS stars. 
^ We derived stellar parameters of the WR stars candidates using a full nonLTE modeling of the observed spectra, 
r""] ' Results. Using a variety of techniques we found that G30 is very young cluster, with age t « 4Myr. The cluster is located in 
I^nJ ' Carina spiral arm, it is deeply embedded in dust and suffers reddening of ~ 10.5±1.1 mag. The distance to the object is 
d=7.2±0.9 kpc. The mass of the cluster members down to 2.35 AIq is ~16OOA40. Cluster's MF for the mass range of 5.6 to 
31.6 A^Q shows a slope of r=— 1.01±0.03. The total mass of the cluster obtained by this MF down to 1 A^0 is about 3 x 10"^ A^q. 
The spectral analysis and the models allow us to conclude that in G30 are at least one Ofpe/WN and two WR stars. The WR 
stars are of WN6-7 hydrogen rich type with progenitor masses more than 60 A^0. 

Conclusions. G30 is a new member of the exquisite family of young Galactic clusters, hosting WR stars. It is a factor of two 
to three less massive than some of the youngest super-massive star clusters like Arches, Quintuplet and Central cluster and is 
their smaller analog. 

Key words. Galaxy: open clusters and associations, stars: Wolf-Rayet, stars: early-type, stars: winds, general-Infrared: general 



1. Introduction 

Stars rarely form in isolation. In fact, it is well known 
that more of the stars in our Galaxy, and in nearby galax- 
ies, are born in groups ranging from small associations 
and open clusters, compact young massive clusters to old 
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globulars. Young clusters are often difficult to find because 
they can be heavily embedded in dust, making them vis- 
ible only in the infrared. We embarked on a long-term 
project to search the inner Milky Way for hidden star clus- 
ters and to study them in details (Borissova et al. 2003, 
2005, 2006; Ivanov et al. 2002, 2005). This project was 
based on the 2MASS (Skrutskie et al. 2006), taking ad- 
vantages of the reduced extinction in the near-IR. Recent 
advances in mid-IR instrumentation made it possible to 
carry out all-sky IR surveys in this spectral region too. The 
recent Spitzer Space Telescope Galactic Legacy Infrared 
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Mid-Plane Survey Extraordinaire (GLIMPSE, Benjamin 
et al. 2003) offers an excellent opportunity to carry out a 
deep census of such objects. GLIMPSE is an excellent tool 
for finding obscured clusters in the Galactic disk because 
the extinction in the mid-IR is a factor of 2-5 lower than in 
the near-IR. Recently, Mercer et al. (2005) implemented 
a comprehensive search for clusters in the mid-IR. They 
used the point source catalog of the GLIMPSE and re- 
ported 92 cluster candidates. In our project we continued 
the investigation of some of these objects using near-IR 
imaging and low resolution IR spectroscopy. 

Massive stars themselves play an important role in 
the ecology of galaxies, providing a major source of ioniz- 
ing UV radiation, mechanical energy and chemical enrich- 
ment. Wolf-Rayet (WR) stars represent an evolved phase 
of the most massive stellar population, and are character- 
ized by high mass loss rates from fast and dense winds. 
Their short lifetimes and high luminosities make them ex- 
cellent tracers of active recent star formation. However, 
serious gaps in our understanding of massive stars exist 
because of the rapid evolution and rarity. WRs in clus- 
ters are particularly interesting because this enables us to 
study a coeval population with uniform metallicity. There 
are about three hundred known WR stars in our Galaxy 
(van der Hucht 2006). At the same time from thousands 
known open clusters and stellar associations small num- 
ber contain WRs and only a few of them host more than 
one such star. This is because a typical cluster (~10'^ A^©) 
contains a limited number of stars more massive than ^^30- 
35 A^0 that can evolve to the WR phase. Therefore, most 
cluster WRs are concentrated in the high-mass clusters, 

1. e. Arches, Quintuplet, Westerlundl. 

Here we report our first results for GLIMPSE 30 (here- 
after G30; Figure[T]) - compact young cluster located near 
the Galactic plane (l = 298?756,6 = -0?408) containing 
many massive stars including at least one Of and three 
Wolf-Rayet members. The presence of at least three WRs 
in G30 puts this cluster into the family of WR rich clus- 
ters. The cluster membership provides excellent observa- 
tional constraints upon the ages and initial masses of this 
type of stars. In this paper we present the main results of 
our photometric and spectral analysis of the cluster and 
the newly discovered WRs. 

2. Observations and Data Reduction 

All observations were obtained with Sofl/NTT (Son of 
ISAAC; Moorwood, Cuby & Lidman 1998). The instru- 
ment is equipped with Hawaii HgCdTe 1024x1024 detec- 
tor, with pixel scale of 0.288 arcsecpx"^. For the spec- 
troscopy we used 1 arcsec slit and the medium-resolution 
grism. The seeing for all observations was 1-1.5 arcsec and 
the sky was photometric. 

Deep JsHKs imaging of G30 were carried out on April 
15, 2006. We took 16 images in each filter in jittering 
mode with 3 arcmin jitter box size to ensure that there 
is no overlapping of the cluster on different images. Each 
individual image was the average of 3x20 sec frames in 




Fig. 1. Pseudo-true color images of the G30. The left 
panel image is composed from Js - blue, H - green and Kg 
- red Sofl/NTT images covering 4f92x4.'92 centered at the 
cluster. The right panel is composed from 3.6 /xm - blue, 
5.8 ^m - green and 8.0 /xm - red Spitzer Space Telescope 
images. The SofI field of view from the left panel is marked 
with a square. The Galactic coordinates I and b are shown 
on this right panel. 

Js, 6x10 sec frames in H, and 10x6 sec frames in Ka,. 
The total integration time was 16min in each filter. To 
obtain photometry of the brightest cluster members we 
took additional shallow images on August 10, 2006, us- 
ing the same jittering pattern but shorter integrations: 
5x1.182 sec. The total integration time was 1.58 min in 
each filter. The data reduction included: sky-subtraction, 
flat-fielding, image alignment and combination into a sin- 
gle final image for each filter. A 3-color composite image 
of G30 is shown on Figure [T] 

The stellar photometry of the final images was carried 
out with ALLSTAR in DAOPHOT II (Stetson 1987). 
The typical photometric errors vary from 0.01 mag for 
stars with Kg ~10mag to 0.10 mag for Ks ^18 mag and 
0.15 mag for Kg ^19 mag. The photometric calibration 
was performed by comparing our instrumental magnitudes 
with the 2 MASS measurements of about 1330 stars, cov- 
ering the color range 0.0< Js— ifs^S.Omag and magnitude 
range 10<Jsrs<14.5mag. A least squares fit of the instru- 
mental jhk magnitudes to the standard 2MASS system 
gave the following relations: 

(Js - Ks) - {j ~k)^ -1.302(±o.oo5) 
Js - j = 0.005(±o.oo7) X (^s - Ks) ~ 2.435(±o.oo3) 
H-h^ -0.021(±o.oo8) X ( Js - Ks) - 2.457(±o.oo3) 
Ks-k^ -0.004(±o.oo6) X (^s - Ks) - 1.134(±o.oo5) 

The final photometry list contains equatorial coordi- 
nates, and JsHKs magnitudes of 7469 stars with photo- 
metric errors less than 0.15 mag. Artificial stars tests show 
that the 80% completeness limit of the photometry is at 
Js = 19.1 and Ks^UA. 

The spectra were obtained on Apr 14, 2006. They cover 
the region from ^2.00 to ~2.35/^m. The slit was aligned 
on a sample of bright cluster members and the telescope 
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Fig. 2. {Js—Ks, Ks) color-magnitude diagrams of 
the stars within the circle with radius R=43.2arcsec 
(R=150px) around G30 (left panel) and a comparison 
field (right panel). See Section [01 for details. 



was nodded along the slit between the exposures. We ob- 
tained 16 images of 150 sec each, on one single slit po- 
sition. In total, spectra of 8 stars were extracted from 
the data. To correct for the telluric absorption we ob- 
served the star HIP59642 (HD106290) - solar near-analog 
of spectral class GIV. The reduction process included: 
sky subtraction, flat flelding, geometric distortion correc- 
tion, image alignment and combination into a final 2- 
dimensional spectrum. Next, we extracted 1-dimensional 
spectra, wavelength calibrated them, and corrected them 
for the telluric absorption by multiplying with the telluric 
standard. After this the target spectra were multiplied by 
a solar spectrum to remove the artificial emission lines due 
to the intrinsic absorption features in the spectra of the 
standard (see Maiolino, Rieke & Rieke 1996). 

3. JsHKs Color-Magnitude and Two-Color 
Diagrams 

3.1. Field-star decontamination 

The cluster G30 is located extremely close to the Galactic 
plane and therefore, in a crowded stellar field. However, it 
is rather concentrated and occupies small area with diam- 
eter ~1.3arcmin, always near the center of the SofI field 
of view. This allows us to define and subtract reliable the 
fore- and background field star contamination. 

Photometric and spatial criteria were used to select 
probable cluster members. First, only the stars in a cir- 
cle with radius R=150px ('^43 arcsec), similar to the 
apparent cluster size, and centered on the cluster were 
chosen. This limits the candidate members to 592 stars. 
The color-magnitude diagram (CMD) of this sample is 
shown on the left panel of Figure[2l The locus of the 
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Fig. 3. {H—Kg, Js—H) two-color diagrams of the stars 
within the circle with radius R=43.2arcsec (R=150px) 
around G30 (left panel) and the comparison field (right 
panel). The continuous and dotted lines represent the se- 
quence of the zero-reddening stars of luminosity classes 
I (Koornneef 1983) and V (Schmidt-Kaler 1982), respec- 
tively. Reddening vectors for 05V and M5I stars are also 
shown. 



cluster Main Sequence (MS) stars can be easily seen at 
1.5<(Js — Ks)<2mag. Some fore- and background stars 
are also present on this diagram. In order to define their 
locus, the CMD of the stars falling into a circle with the 
same radius but centered lOOarcsec to the North-East 
from the cluster's center is shown on the right panel of 
FigureO The contamination become significant only at 
K > 17 mag while the rest of the cluster locus is almost 
entirely devoted of field stars. 

The majority of the cluster candidate members oc- 
cupy a well-defined locus on the two-color diagram (TCD) 
as well, at H—Ks^O.Gmag and Js — H^l.l mag. Figure[3] 
shows that there is only small number of field stars falling 
in the cluster locus on the right panel of the same figure 
where the TCD of the comparison field (defined above) is 
shown. 

The decontamination procedure was based on the color 
and magnitude distributions of the field stars in five stel- 
lar fields around G30, to improve the field star statistics. 
Each field covers the same area as the cluster. Their com- 
bined CMD was divided into bins and the number of stars 
in each bin was divided by five to normalize the combined 
field area to the cluster area. Then, we randomly sub- 
tracted from the corresponding bin of the cluster's CMD 
as many stars as were present in the corresponding field 
CMD bin. The bin sizes varied depending on the number 
of stars per bin. In other words, we merged nearby bins, if 
they contained less than 2 stars. The total number of stars 
subtracted throughout the decontamination was 191. 
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Table 1. Of and WR stars in G30. 



ID 


02000 


^2000 


Js-H 






1 


12:14:33.91 


-62:58:48.7 


1.250 


0.732 


9.676 


2 


12:14:33.09 


-62:58:51.0 


1.194 


0.648 


9.674 


3 


12:14:31.76 


-62:58:51.9 


1.188 


0.577 


8.625 


4 


12:14:31.54 


-62:58:54.3 


1.152 


0.573 


8.605 



3.2. Analysis of the decontaminated diagrams 

The decontaminated CMD and TCD of the cluster are 
shown in FigureH) The CMD morphology presents an ex- 
tended MS reaching down to ^^ITmag. The majority of 
these stars form well-defined compact sequence on the 
TCD around H—Ks^O.Gmag. The spectroscopically con- 
firmed Of/WN and WR stars (marked with diamonds) are 
among the brightest stars in the field and they are located 
to the right of MS, as expected for evolved objects. There 
is a hint for deviation of the most massive stars from the 
MS at 10<Xs<12mag. 

The CMD shows a population of 45 likely pre-main 
sequence stars (PMS; marked with crossed circles) at 
l<ii's<17mag and 2< Js — Js<2.5 mag. We can not ex- 
clude the possibility that some of them might be unre- 
solved binaries or MS stars subjected to higher differen- 
tial reddening. To constrain the age and masses of these 
stars we applied isochrone fit with PMS tracks of differ- 
ent ages: 0.1, 1, 4, 7, and lOMyr from Siess, Dufour & 
Forestini (2000) shifted with the obtained distance modu- 
lus and reddening to the cluster (see Section^] for details). 
Considering the photometric errors and the crowding ef- 
fects we reach tentative conclusion that the majority of 
the PMS candidates have ages between 1 and 10 Myr. 

Non-negligible number of objects (23) fall to the right 
of the reddening vector for the reddest stars on TCD. 
Most likely they have excess due to circumstellar envelopes 
and/or discs. We marked with sohd dots those stars that 
are more than 3cr ^0.12 mag redder than the reddening 
vector for the reddest MS stars. 



4. Wolf-Rayet Stars 

The Ofpe/WNL and WR members of G30 are marked 
with diamonds in Figure |4]). Their coordinates and in- 
frared broadband properties are listed in Table[l] and their 
spectra are plotted in Figure[5l The S/N at A=2.07/im for 
the stars 1 to 4 is 90, 80, 80, and 60, respectively. 

The equivalent widths of the emission lines were mea- 
sured from the normalized spectra. To quantify the sub- 
types of WR stars, we compared our data with the near- 
IR spectral atlases of Crowther & Smith (1996) and Figer, 
McLean & Najaro (1997), following the spectral classifi- 
cation scheme of Crowther et al. (2006). The results are 
summarized in Table [2l 
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Fig. 5. Infrared spectra of the G30 stars. The spectra are 
continuum normalized and shifted vertically by 0.5 for dis- 
play purposes. The S/N at A=2.07/im for the stars 1 to 4 
is approximately 90, 80, 80, and 60, respectively. 

Our analysis shows that star No. 1 is WN7, based on 
the presence of strong N III at 2.116 fim. The He II line at 
2.189 /im is weaker than Br7, but He I lines are also weak. 

Star No. 2 is WN6 because the observed ratio He II 
2.189 /im / Br7 > 1 shows lack of hydrogen. The spectrum 
is similar to that of WR136 = HD192163 (van der Hucht 
2001). 

Star No. 3 is particularly interesting. The presence of 
Br7 emission, the weak He I lines and the lack of Hell 
lines neither in absorption nor in emission suggest that it 
is Ofpe/WNL. The spectrum resembles these of HD 47129 
(08p), HD 148937 (06fp), HD 152408 (08Iafpe), and 
HD 269582 (0fpe/WN9). There are only several known 
stars of this type in the Galaxy. Interestingly, Ofpe stars 
also reside in clusters: in Quintuplet (Figer, McLean & 
Morris 1999) and in the Galactic Center cluster (Cotera 
et al. 1999). This warrants further study of the star. 

Star No. 4 shows strong Hell line (as strong as Br7), 
but it also has He I and no N III line. The measured line 
intensities (Table [2|) suggest a tentative WN7-8 classifica- 
tion. Unfortunately, the star was located at the edge of the 
slit and its spectrum has the lowest S/N ratio of all four 
WR and O star candidates. We compared this spectrum 
with the template spectra of WN and O stars given in 
Conti et al. (1995), Hanson, Conti & Rieke (1996), Figer, 
McLean & Najaro (1997) and Crowther et al. (2006), 
smoothing the templates to the resolution of our data. The 
closest matches are with WR131 (WN-Ha) and HD16691 
(04If+), shown in Figure[Sl Note that the intensities of 
He II and Br7 are higher in the template spectra and they 




Jg-Kg H-Kg 

Fig. 4. Left: Decontaminated {J~Ks, Ks) color-magnitude diagram. An isochrone for 4 Myr from Girardi et al. (2002) 
and PMS isochrones for various ages and solar metallicity from (Siess, Dufour & Forestini 2000) are shown. The 80% 
completeness limit of the photometry is marked with a dashed line. Solid dots indicate stars with color excess, as 
defined on the right panel. Right: Decontaminated {H—Ks, Js—H) two-color diagram. The continuous and dotted 
lines represent the sequence of the zero-reddening stars of luminosity classes I (Koornneef 1983) and V (Schmidt-Kaler 
1982), respectively. Reddening vectors for 05V and M5I stars are also shown. The probable PMS stars are marked 
with crosses. The 3cr-limit line of the reddening vector for OV star is drawn «0.12mag to the right of the reddening 
vector. The stars redder than this line show color excess and are marked with solid circles. 



have much more reliable Hcl + NIII line. This compar- 
ison also shows similarities with WN7 class but the Hell 
emission suggest early or mid-0 star. 

Another four stars fall into the slit by serendipitously 
but their S /N is too low for quantitative analysis and they 
show no strong emission lines. 

This spectral classification described above gives a 
rough estimate of the temperature and luminosity of the 
stars. The precision of these estimates is related to the 
intrinsic variation within the classes. Unfortunately, the 
parameters of the WR stars are not truly homogeneous. 
Later in the text we use the WR stars as distance indica- 
tors, so having better estimates is crucial. 

The limited number of emission lines in the K atmo- 
spheric window makes the line ratios less reliable than a 
full line modeling. We used the CMFGEN model (Hillier & 
Miller 1998) to obtain synthetic line profiles. Even though 
the spectra do not contain a lot of spectral features we 



run models with complex atoms including CNO, Si, P, S 
and Fe because they have effect on the electron density, 
temperature, etc. and therefore - on the properties of the 
Hydrogen and Helium lines. 

The two main lines in our spectra Br7 and 
He II 2. 18 /im are sensible to both temperature and mass 
loss rate. Their ratio is sensitive to the helium abundance. 
To measure all parameters simultaneously we calculated 
a grid of models spanning a range of temperatures from 
35000 to 40000 K and mass loss rates from 1x10"'^ to 
1.5x10^^^ A^/yr. The chemical composition was set to so- 
lar except for He/H which was set to He/H=1.0, in terms 
of number of atoms. Initially, the luminosity of every star 
was set to logL/LQ=6.0 and the terminal velocity to 
Voo = 1000km/s, determined from the width of the spec- 
tral lines. There are no lines with good P Cyg profile in 
the K spectrum, so the value of Voo should be treated with 
caution. We discus the values of the luminosity below. 
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Table 2. Spectral classification of Of and WR stars in G30. For every object the equivalent widths (first hne) and 
FWHM (second line) in A for prominent near-IR lines in the spectrum are presented. The central wavelengths of the 
lines are in fim. 



Star 


He I 


NV 


Hel + NHI 


He n + Br7 


HeH 


HeII/Br7 


Type 




2.058 


2.110 


2.115 


2.165 


2.189 






1 


PCyg 


no 


8.97 


23.34 


14.84 


0.67 


WN7 








97.2 


68.7 


60.7 




weak 


2 


PCyg 


no 


-2.6 


10.56 


15.36 


1.45 


WN6 








82 


74 


89 




weak 


3 


PCyg 


0.39 


0.98 


6.26 


no 


0.67 


Ofp/WN 






19 


30 


47 






? 


4 


absorb. 


no 


no 


2.77 


2.03 


0.83 


WN7-8 




-2.53 






45 


49 




weak 



1.4 























\ A HD'6691 (0-1 If*) 





2.0 2.1 2.2 2.3 2.4 

Wavelength (nm) 

Fig. 6. Comparison of the spectrum of star #4 with tem- 
plate spectra of WR131 (WN-|-a) and HD16691 (04If+) 
from Figer, McLean & Najaro (1997) and Hanson, Conti 
& Rieke (1996). 



The synthetic spectra were degraded to the resolu- 
tion of the observed spectra, using the profile of the arc 
lamps as a measure of the instrumental profile (Figure[7| . 
Then, we measured the intensities of Br7, He 1 2.05 /im, 
HeI2.il/im and He II 2.18 /im on both the models and 
the observed spectra. The grid of models defines a sur- 
face in Teff vs. M space. The observed value of every 
spectral feature is a iso-line on that surface. The crossing 
of the iso-lines, determined by different features, sets the 
parameters which reproduce simultaneously all measured 
features. The temperature depends mainly on the ratio of 
He I/He II lines while the intensity of Br7 is mainly related 
to M. Finally, the intensity of He II 2.18 iiui is sensitive to 
Teff, Ai and to the helium composition but the first two 
are already constrained by the other parameters. We had 
to adjust the He/H ratio to reproduce this line well. 




2.00 2.05 2.10 2.15 2.20 2.25 

Wavelength (|im) 

Fig. 7. CMFGEN model (Hilher & Miller 1998) fit to the 
observed spectra of stars #1, #2, and #4. 



As pointed out above, the values of the stellar lumi- 
nosities were assumed. There is no easy way to determine 
the exact luminosity from the spectra alone. We broke the 
luminosity-radius-mass loss rate degeneracy applying the 
relation between modified wind momentum H=A^VooR2'^ 
and the luminosity (Kudritzki, Lenon & Puis 1995). The 
derived Teff, M and Voo set a family of possible mod- 
ified wind momenta and stellar luminosities, where the 
stellar radius is a free parameter. We calculated the pos- 
sible n and L for a range of radii and plotted them on 
the Kudritzki, Lenon & Puis (1995) diagram together 
with data for WNha stars from Hamann, Griifener & 
Liermann (2006). The crossing between this line and the 
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Fig. 8. Modified wind momentum 11 versus luminosity L 
diagram for a range of radii, for stars #1, #2, and #4. 
The crosspoint between the individual Hne for every star 
and the observed log L— log 11 relation determines the stel- 
lar luminosity. WNha stars from Hamann, Grafener & 
Liermann (2006) are marked with diamonds. The least 
square fit to these data alone is shown with a dashed line. 
O stars from Lamers et al. (1999) are marked with squares. 
The least square fit to all data is shown with a solid line. 



observed log L— log 11 relation determines the stellar lumi- 
nosity which gives the same emitted spectrum and satisfies 
the log L— log n relation - the least square fit to the data 
(open diamonds) and is shown in Figure[8] with a dashed 
line. The obtained logL/L© is 5.77, 5.57 and 5.43 for stars 
#1, #2, and #4 respectively. However, the sample of WN 
stars, containing hydrogen is very small and we added the 
data of 0-type stars given in Lamers et al. (1999; shown 
with squares) to improve the statistics. The new relation 
based on the expanded dataset is shown with continuous 
line in Figure[51 The luminosities are listed in Table[3] and 
we will use them as our final estimates. 

This technique allowed us to determine the tempera- 
ture and the mass loss rate of the strong Hell2.18/xm 
emission stars #1, #2, and #4. We cross-checked the 
model predictions with the features expected in K spec- 
tra. The Nlll at 2.115/im is heavily blended with a Hel 
line which makes the accurate determination of the nitro- 
gen composition difiicult. Nevertheless, the profile fitting 
favors an increase of N abundance by a factor of 10 with re- 
spect to the solar. Higher resolution spectra are required 
to achieve better precision. The CIV 2.07 /zm is absent 
from the observed spectra, which suggests a reduction of 
the carbon abundance by a factor of at least 10. In gen- 
eral, these three stars have increased mass loss rate and 
relatively low terminal velocity, higher helium abundance 
and reduced carbon. All these characteristics make the 
first two of them good candidates for hydrogen rich WN 
stars. The type of the third star can note be constrained 
with our data. 



Table 3. Physical parameters of WR stars in G30. The 
bolometric corrections BC from Crowther et al. (2006), 
the absolute K magnitude and the true distance modulus, 
obtained from WR stars are also listed. 
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logL 




log 


H/He 


BC 


Mk 


{K-Mk)o 




\Lo\ 


[K] 


[Me/yr] 




[mag] 


[mag] 




1 


6.04 


37700 


-4.91 


1.000 


-3.8 


-6.55 


15.065 


2 


5.91 


38185 


-5.12 


1.030 


-4.4 


-5.62 


14.134 


4 


5.81 


36800 


-5.30 


1.015 


-3.6 


-6.18 


13.625 



5. Discussion 

5.1. Reddening and Distance 

We used photometric and positional criteria to se- 
lect probable cluster members. The candidates oc- 
cupy a well-defined locus on the color-color diagram 
at H-Ks~Q.b-0.& mag and Js--ff~0.9-1.2 mag 
allowing us to determine the reddening to the clus- 
ter. First, we simply measured the color excesses of 
this locus on the color-color diagram with respect 
to the sequences of unreddened MS stars (Schmidt- 
Kaler 1982), obtaining E(iJ-ii's)=0.60±0.06mag 
and E(Js— i?) = 1.10±0.06mag, corresponding to 
Ai<-s=1.12±0.08mag and Ay=10.00±0.7mag. Through 
the paper we used the reddening law of Rieke & Lebofsky 
(1985). 

The reddening can be determined also from the 
apparent MS colors, given that the MS of young 
clusters is nearly vertical. Applying this method and 
fitting only unevolved cluster MS stars with the 
4Myr Geneva isochrone (Girardi et al. 2002) we 
obtained E(Js— i?) = 1.81±0.12mag, corresponding to 
A/<-s=1.19±0.12mag and Av = 10.65±l.l mag. 

We verified these estimations once again against 
the spectral type classification of the WR clus- 
ter members, adopting the intrinsic WR colors of 
Crowther et al. (2006). The averaged color excesses 
for the WN stars are E(iJ-ii:s)=0.56±0.07mag 
and E(Js— i^s)=l-79±0.13mag, corresponding to 
Aifg=1.18±0.13mag and y4.y=10.54±1.3mag, in reason- 
able agreement with our previous estimates. 

Ideally, to obtain an accurate distance to the clus- 
ter we need the spectral class and respectively the ab- 
solute magnitude of some of the MS stars because they 
are much more uniform in comparison with the WR stars. 
Unfortunately, we do not have spectra of MS stars. The 
WN luminosities however allow us to obtain relatively 
good estimate of the distance. The bolometric corrections 
(BC) for WN stars from Crowther et al. (2006) were used 
to transform the luminosities to absolute Ka^ band magni- 
tudes. The BCs and the distance moduli of G30 WR stars 
are listed in Table[3] Averaging the individual estimates, 
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we obtain A/fs=1.16 and (Xg— Afxs)o = 14.31±0.35mag, 
corresponding to distance of d=7.2±0.9kpc. 

We also used the mean absolute iiTg-band magnitudes 
of the WN subtypes given in Crowther et al. (2006) to ver- 
ify our distance modulus estimate. Averaging over the G30 
WR members we obtained (i^s— Afi<'s)o=13.58±0.35mag. 
This corresponds to a shorter distance of d=5.2±0.9kpc. 

In our further analysis we will use the larger dis- 
tance modulus to the cluster because the models of WR 
stars yield more reliable temperature and luminosity es- 
timates (and respectively - distance) than the spectral 
classification based on comparison with template spectra. 
Furthermore, the shorter distance modulus moves the hy- 
drogen burning turn-on point, (where the isochrone of the 
PMS of 4Myr reaches the MS), to about 1 mag above 
the observed lower MS cut-off. This is in disagreement 
with the estimated age of the cluster (see Figured and 
Subsection IS . 21 for details). 

5.2. Cluster Age 

It is difficult to obtain a reliable age using only an 
isochrone MS fitting of young stellar clusters because of 
the nearly vertical linear MS locus. The first constrain 
comes from the lack of red supergiants in G30, evident 
from the CMD (Figure [2]). The stellar evolutionary mod- 
els with rotation predict the onset of red supergiants at 
~4.5— 5Myr, defining an upper limit to the cluster age. 

The WR phase is very short lived and the presence of 
WR stars limits the maximal age of the cluster (Meynet 
& Maeder 2005). AU known WR stars in G30 are of the 
WN6-7ha (hydrogen rich) subtype. There are indications 
that hydrogen rich WN7 stars are descendants of massive 
stars with initial masses above 50— 607Wo (e.g., Crowther 
et al. 1995). Then an upper age limit of 4 — 4.5Myr can be 
set, independently of the exact metallicity and mass- loss 
scenario. 

The MS turn-off MS point provides a con- 
sistency check. Our final true distance modulus 
(i<rs— -Mi<'s)o=14.31mag suggests that the brightest un- 
evolved star is 09.5 or BO setting the cluster age to 
~3— 4.5Myr. This estimate is not independent, because 
we used WR stars to obtain the distance to the cluster. 

The age of the PMS stars was determined by fitting 
theoretical 0.1, 1.0, 4.0, 7.0, and lOMyr PMS isochrones 
from Siess, Dufour & Forestini (2000) to the CMD 
(Figured]). Note that the PMS stars spread over wide 
age range but the main locus is between the 1 and 4 Myr 
isochrones. There are some stars with ages less than 1 Myr 
and also a possible concentration of stars around 10 Myr 
but it is too close to the depth of our photometry to draw 
a certain conclusion. We refrain from making more accu- 
rate conclusion about the age and the age spread because 
a significant fraction of the PMS's falls into the zone of 
photometric incompleteness. Deeper CMD is needed for 
that. However, we point that continuous star formation 



scenario or at least an extended burst can not be excluded 
to have occurred in G30, based on our data. 

The fraction of stars with NIR-excess correlates in- 
versely with the stellar age, over small age ranges 
(Hillenbrand 2005). The vast majority (~90%) of stars 
older than 3-8 Myr ceases to show evidence for accretion. 
The fraction of IR-excess stars in a very young stellar clus- 
ter such as G30 can be used as an age indicator. We deter- 
mined that for G30 the fraction is 12%. The empirical cal- 
ibration of Hillenbrand (2005) suggests an age of 3-4 Myr, 
in perfect agreement with our previous estimations. 

Similar to Hillenbrand, Bauermeister & While (2007) 
we point out that ages may be affected by the photometric 
uncertainties and astrophysical effects such as variability 
of young objects, unresolved binaries, etc. In this regard 
it is important to discuss the following two questions: 

1) How the distance uncertainty affects the PMS-age? 
- As discussed above, we have two distance moduli esti- 
mates that differ by w 0.7 mag. Propagated to the PMS 
age this corresponds to age difference of 6-7 Myr. The 
shorter distance modulus gives an age of more than 10 Myr 
which is inconsistent with the presence of WNha stars so 
the true age uncertainty is even smaller. 

2) How the photometric errors affect the age spread? - 
The formal photometric uncertainties at the turn-on point 
are ct( Jg— i^s)=0.12 and CT(i4rs)=0.10 mag. Adding the er- 
rors due to the transformations to the standard JgHKs 
system, we see that the contribution of the photometric 
errors to the age range does not exceed 1-2 Myr. Therefore, 
most of the observed age spread is intrinsic to the PMS 
population. 

5.3. Initial Mass Function and Total Cluster Mass 

To calculate the IMF we converted the stellar magni- 
tudes into masses using the Geneva models. The to- 
tal mass of the observed cluster members down to 
--2.35 A^0 (including the WR stars) is ^1600 Mq adopt- 
ing a distance modulus (iiTs— -^i<:s)o = 14.31 and redden- 
ing Aifs=1.16mag. The mass estimation of G30 mea- 
sured from the {{Js—Ks)o, M^g) diagram as described 
in Borissova et al. (2003; see their Fig. 8) gives similar 
result. 

We also constructed a background subtracted mass 
function of the cluster and fitted it with a sin- 
gle power-law, obtaining F^— 1.01±0.03 (in this terms 
the Salpeter slope is r = — 1.35) over the mass range 
logTW/TWo = 0.75 - 1.5 (5.62-31.627^0). The integra- 
tion over this MF down to I AAq combined with the 
masses of the WR stars leads to a total cluster mass of 
~3 X 10'^ A^0. Naturally, this is only a lower limit. 

These results make G30 only two to three times less 
massive than some of the most massive young clusters in 
the Galaxy - Arches and Quintuplet. 
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5.4. Comparison with Other Massive Clusters and 
Position in the Galaxy 

The physical properties of the G30 together with data 
for Quintuplet, Arches and the Central cluster from Figer 
(2004) are summarized in TablejH Here A^i is the total 
cluster mass in observed stars, A^2 is the total cluster 
mass for all stars integrated down to IMq, assuming a 
Salpeter IMF. The "Radius" is the radius in parsecs, from 
the central stellar surface density peak. The mass densities 
pi and p2 are simply A4i and A^2 divided by the cluster 
volume. The "Age" and "Luminosity" are the estimated 
cluster age and total luminosity. 

Quintuplet shows very similar characteristics to G30. 
They both have bigger radius and lower stellar den- 
sity in comparison with Arches and Central cluster but 
Quintuplet is a factor two more massive and a factor of 
three more luminous than G30. The Quintuplet is also 
older than G30 because it contains one RSG. The proper- 
ties of G30, especially the presence of WR stars, make G30 
a smaller analog of Arches, Westerlund 1 and Quintuplet. 

The galactocentric distance of G30 is i?GC=8.1kpc, 
assuming that the Sun is at i?GC, 0=8-4 kpc from the 
Galactic Center (Figure[9|) . G30 is located between Carina 
and Crux spiral arms, closer to Carina. It probably be- 
longs to some lateral branch of the farther inner side of 
this arm. We see the cluster through a dust window in the 
nearer side of the Carina arm. Note that if the distance to 
the cluster is underestimated it can belong to the farther 
side of the main Carina arm. 



-180° -150° 180° 150° 120° 

\^ ^ \ \ \ I / / / 7 7 
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Fig. 9. Location of G30 in the Galaxy. The galactocentric 
distance of the cluster is Rgc=8.1±0.4 kpc. The error 
bar shows the la error of 0.9 kpc in the distance to the 
Sun. The spiral structure is taken from Cordes & Lazio 
(2002). The solid lines show the modified spiral model of 
the Galaxy from Taylor & Cordes (1993). The dashed lines 
are a four-arm logarithmic spiral model combined with a 
local arm. The "-I-" sign marks the Galactic center and 
the Sun is labeled. 



6. Summary 

We report new results of our long term project to study ob- 
scured Milky Way star clusters. We obtained deep JsHKs 
imaging of G30, a dense stellar cluster with spectroscop- 
ically confirmed WR stars and a sizable population of 
young stellar objects. G30 is deeply embedded into gas and 
dust and suffers a reddening of Ay ~ 10.5 mag. The object 
probably belongs to the Carina spiral arm and is located 
at a distance of 7.2±0.9kpc from the Sun. The cluster is 
approximately 4Myr old. The uppermost MS stars have 
evolved away from the zero-age MS. G30 is massive, with 
a lower limit of the total mass of 3 10^ Mq. This estimate 
includes only stars with masses above ~1.OA^0. 

Spectral analysis and modeling of K spectra for four 
objects show that one of these is Ofpe/WN star, two are 
hydrogen rich WN6-7 stars, and the last is a WN or 0-type 
star, all with progenitor masses above 60 A4q. The CMD 
suggests that there might be more WR or O type cluster 
members and additional observation are planned address 
this possibility. G30 is a new member of the exquisite fam- 
ily of massive young Galactic clusters, hosting WR stars. 
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Table 4. Parameters of G30 compared to the massive clusters in the Galactic Center (see Section [5751 for more details. 
The data for Quintuplet, Arches and the Central cluster are from Figer (2004; see their Table 1). Note that the total 
mass and the corresponding density of G30 are only lower limits. 
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